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ABSTRACT 

Blazars exhibit flares across the electromagnetic spectrum. Many 7 -ray flares are highly correlated 
with flares detected at optical wavelengths; however, a small subset appears to occur in isolation, with 
little or no variability detected at longer wavelengths. These “orphan” 7 -ray flares challenge current 
models of blazar variability, most of which are unable to reproduce this type of behavior . We present 
numerical calculations of the time-variable emission of a blazar based on a proposal by |Marsclier et ah] 
to explain such events. In this model, a plasmoid (“blob”) propagates relativistically along the spine 
of a blazar jet and passes through a synchrotron-emitting ring of electrons representing a shocked 
portion of the jet sheath. This ring supplies a source of seed photons that are inverse-Compton 
scattered by the electrons in the moving blob. The model includes the effects of radiative cooling, a 
spatially varying magnetic field, and acceleration of the blob’s bulk velocity. Synthetic light curves 
produced by our model are compared to the observed light curves from an orphan flare that was 
coincident with the passage of a superluminal knot through the inner jet of the blazar PKS 1510—089. 

In addition, we present Very Long Baseline Array polarimetric observations that point to the existence 
of a jet sheath in PKS 1510—089, thus providing further observational support for the plausibility of 
our model. An estimate of the bolometric luminosity of the sheath within PKS 1510—089 is made, 
yielding L s h ~ 3 x 10 45 erg s _1 . This indicates that the sheath within PKS 1510—089 is potentially 
a very important source of seed photons. 

Subject headings: galaxies: active - galaxies: jets - polarization radiation mechanisms: non-thermal 
- relativistic processes - techniques: interferometric 


1. INTRODUCTION 
Blazars are the most luminous pe rsiste nt objects in the 


sky (up to ~ 10 48 erg s 4 ; e.g., Krolik 1999). They ex¬ 


hibit variability across the electromagnetic spectru m on 


timescales ranging from minutes to years (e.g., Aharo- 
nian et al. 2007 and references therein). Blazars consti- 
tute a sub-class of active galactic nuclei (AGNs) whose 
relativistic plasma jet s are t hought to be closely aligned 
to our line of sight (Blandford & Konigl 1979). The 
emission from blazars is predominantly non-thermal, in¬ 
cluding radio synchrotron radiation and X-ray and 7 - 
ray inverse-Compton scattering. The central engines of 
blazars, which are believed to be the ultimate source 
of this variability, are unresolved and opaque to radio 
waves. Hence, they are not imaged directly, even with 
the 0.1 mas resolution obtained with the Very Long Base¬ 
line Array (VLBA) at 43 GHz. Despite this limitation, 
by monitoring and subsequently modeling flares in the 
high-energy emission from these AGN, we can potentially 
gain insight into the sub-parsec-scale physics of the jets 
close to the central engines. 

Many theoretical models have been developed to ex¬ 
plain both the spectral energy distributions (SEDs) and 
high-energy variability of blazars. The SEDs of most 
blazars take the form of two humps - one in the ra- 
dio/UV portion of the spectrum (synchrotron) and the 
other at higher energies in the X-ray/ 7 -ray portion of 
the spe ctrum (interpreted her e as inv e rse-Compton emis- 


sion 


Blandford fc Konig l (11979) 


(1985), and Hughes, Ailer, & Ailer (1 


Marscher & Gear 
proposed that 


shock acceleration within the jet could result in the ob¬ 
served varia tions in radio emission. IDermer fe Schlick-I 
(1993) presented a leptonic model in which the 


eiser 


high-energy hump of the SED is produced by comp- 
tonization of accretion disk photons by electrons con¬ 
tained within a plasmoid (“ blob”) moving relativ istically 
along the jet. Alternatively, Bottclier et al. (2009) formu¬ 
lated a hadronic model in which the high-energy hump 
is produced through py pion production from relativis¬ 
tic protons within the jet encou ntering bot h in tern al 
and ex ternal photon fields. |Sikora, Begelman, fc Rees) 
(1994| also explored the effects of comptonization of ex- 
ternaf radiation fields from the accretion disk and the 
surrounding broad e missio n-line region (BLR) by elec¬ 


trons within the jet. Spada et al. (2001) presented an 


internal shock model for blazar variability in which shells 
of plasma within the jet collide, producing forward and 
reverse shocks that result in particle acceleration within 
the jet. While these models can explain many aspects of 
the radiative mechanisms potentially at work in blazars, 
there remain aspects of high-energy variability that are 
not accounted for by these theories. 

Recent observations obtained with the Fermi Large 
Area Telescope (LAT) and the VLB A have shown that a 
large number of 7 -ray flares in blazars are strongly corre¬ 
lated with the passage of sup erluminal knots throu gh the 
radio cores of these objects (Marscher et al.||2012|. The 
majority of these 7 -ray flares tend to be associated with 
flares detected at optical/UV wavelengths. The high ca¬ 
dence of observations provided by the Fermi spacecraft, 
however, has also allowed the identification of a small 
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Fig. 1. — 7 -ray (upper panel) and optical (lower panel) light 
curves of PKS 1510—089 during 2009. The dashed vertical lines 
delineate an example of an orphan 7 -ray flare in which the promi¬ 
nent flare in the GeV band has at most a weak counterpart in the 
optical band. The vertical arrow in the lower panel marks the time 
when a superluminal knot passed through the 43 GHz core of PKS 
1510—089, coincident with a large optical/ 7 -ray fla re that occurred 
rough ly 20 days after the onset of the orphan flare ||Marscher et al.| 
[2^010). The horizontal bar represents the uncertainty m tfuTtime 
oi the knot’s passage through the radio core. 


sub-set of 7 -flares that seem to occur in relative isolation, 
with little or no variability detected in the other bands. 
These isolated events are termed “orphan” flares. In one 


case, 


Marscher et al. (2010) present observations of an 


orphan 7 -ray llare that occurred during the passage of 
a superluminal knot through the inner jet of the blazar 
PKS 1510-089 (see Figure jll 
Based on the orphan 7 -raynare highlighted in Figure[l] 
here we develop a model of blazar variability to account 
for this type of behavior. In our model, a blob consist¬ 
ing of a power-law distribution of electrons propagates 
relativistically along the jet axis of the blazar and passes 
through a synchrotron-emitting ring of electrons, which 
represents a shocked segment of a more slowly moving 
jet sheath (see Figure [2] and the sketch presented in Fig¬ 
ured] below). The ring creates a very localized source of 
seeaphotons that are inverse-Compton scattered by the 
electrons in the moving blob, creating an orphan 7 -ray 
flare as the blob passes through the ring. We point out 


(see e.g., 

Ghisellini & Madau 1996, 

Bottcher & Dermer 

1998 and 

B6ttcher||2005p that require a reflecting cloud 


close to the jet axis ro produce 7 -ray flares. 

The presence of a non- to mildly-relativistic sheath of 
plasma enshrouding the relativistic spine of a blazar jet 
has been discu ssed by sever al auth o rs. From a theoreti- 
cal standpoint, Ghisellini, Taveccliio, & Chiaberge (2005) 



Fig. 2. — Proposed model for an orphan 7 -ray flare: a blob 
containing relativistic electrons with a power-law distribution of 
energies passes through a synchrotron emitting ring and inverse- 
Compton scatters the ring photons. 


investigated the radiative interplay between these two re¬ 
gions within the jet via the inverse-Compton mechanism. 
More recently, the SEDs of several blazars have been fit 
wit h models that utilize th e presence of a jet sheath (see , 
e.g., Aleksic e t al.|2014 and Thvecch^^ 


In addition, Tavecchio, Ghisellini, & Guettal2UI4) dis- 


cuss the potential of sheaths of plasma within blazars 
to produce the high-energy neutrinos being detected by 
IceCube. Several authors have also discussed what the 
observational signature o f a sheath of plasma wo uld be 


within a blazar jet (see Attridge, Roberts, & Wardle 


1999). Velocity shear between the jet and the ambient 


medium into which the jet propagates is expected to align 
the magnetic held on the outer edges of the jet with the 
jet axis so that the sheath polarization should be high 
(tens of percent) relative to the spine. The electric vector 
position angles (EVPAs) associated with this polariza¬ 
tion should be perpendicular to the jet axis. This polari- 
metric signatu re has indeed been dete cted in a number 


of blazars (see Pushkarev et al. 
larimetric observations of PKS 


20051. We present po- 
510—089 that point to 


the existence of a jet sheath within this blazar, lending 
support to our model of orphan 7 -ray flares. 

This paper is organized as follows: In §2 we outline 
the theory relating to the radiative transfer modeled in 
our code. In §3 we describe the numerical simulation 
performed and list all of our model parameters. In §4 we 
compare synthetic light curves produced by the model to 
the observed optical and 7 -ray light curves of the orphan 
flare of PKS 1510—089. In §5 we present VLBA polari- 
rnetric observations that point to the existence of a jet 
sheath within PKS 1510—089. In §6 we present our con¬ 
clusions and summary. We adopt the following cosmo¬ 
logical parameters: H 0 = 70 km s -1 Mpc -1 , D m = 0.3, 
and 11 a = 0.7. 


2. THEORY 

Models of blazar emission, whether leptonic or 
hadronic, rely on the presence of an emission region (a 
“blob”) moving relativistically along the jet. For the 
purposes of this work we adopt the leptonic scenario and 
treat the blob (and ring) as consisting of electrons with 
a power-law energy distribution, 71,2(7) oc j~ s , over the 




















































7 -ray Variability In Blazars 


3 


range 7 m i n to 7 max , where the energy of a given electron 
is E = -fm e c 2 . 

There are three main components of the non-thermal 
emission produced by the passage of a blob of plasma 
down the spine of a jet in our model. The first is syn¬ 
chrotron radiation (syn) produced by electrons gyrat¬ 
ing around magnetic field lines within the blob. These 
electrons also inverse-Compton scatter the synchrotron 
photons they produce, a process known as synchrotron 
self-Compton (ssc) emission. The third component is 
inverse-Compton scattering of the external (to the blob) 
ring synchrotron photons by the blob’s electrons, a pro¬ 
cess known as external Compton scattering (ec). As the 
blob approaches the ring, the photon density in the co¬ 
moving frame of the plasma increases, resulting in an 
ec driven 7-ray flare that peaks and decays as the blob 
passes through and then moves away from the ring. 

The photon production rates (photons s -1 cm~ 3 e _1 ) 
for these three physical processes, n syn , n ssc , and fi ec , 
are given below (e = hv/m e c 2 is defined as a dimension¬ 
less photon energy). In our radiative transfer code, we 
evaluate these photon production rates in the co-moving 
frame of the blob, after which a series of Lorentz trans¬ 
formations is applied to obtain the resultant flux in the 
observer’s frame. In all of the following equations, the 
dimensionless electron velocity, /3 e = v e /c ss 1, has been 
set to unity. 

The synchrotron photon production rate per unit vol¬ 
ume for electrons within the blob or the ring is given 
by: 


<, syn 


(e) = 


y/3e 3 B 
27 jhm P c 2 


F(x) n e ( 7) d7 (1) 


(e.g., Joshi & Bottcher 20111, where x = v / v c is the ratio 
of the observed frequency to the critical frequency (77), 
which is given by 17 = 3ej?7 2 /47rm P c ; and F(x ) is an 


appro ximation to a Bessel function (see Joshi & Bottcher 


2011 ). 

The ssc photon production rate per unit volume within 
the blob is given by: 


M = 


47r 


3(7) d 7 


« P h° ( e ) S(e s ,e,7) de 


( 2 ) 

(e.g., Joshi & Bottcher 2011), where n b [“ b (e) = 

h'syn(e) fbiob/c is the synchrotron photon density pro¬ 
duced by the electrons in the co-moving frame of the 
blob. The isotropic scattering cross-section, g(e s , 6,7), is 
given by: 


2(e s ,e,7) = 


3co~t 

167 4 e 


(^- 1 ) 


47 2 


4e7 


[Cl + C2 + C3] for e < e s < 1+4e7 > 


where a t is the Thomson cross-section, e is the incoming 
photon energy, e s is the scattered photon energy, and Ci, 
C2, and C3 are defined as follows: 


Ci = 277 In (77) 

C2 = (1 T 277) (1 — 77) 

_ (4e7?7) 2 (1 - 77) 
^ 3_ 2(l+4e 7 7 ? ) 


(3) 


where 77 = e s /4e7 2 (l — e s /j) ( Joshi fe B6ttcherj|201l| ). 

The general ec photon production rate per unit volume 
within the blob for the scenario presented in Figure [2] is 
given by: 


/ /* e max (‘ /*7max 

d ^ / 7ip“ s (e) de ®dfl e / n e { 7) 

6m i n J J 7m i n 


x (1 — cos^) d7 


dcr 


dCad^c 


( 4 ) 


(Dernier & Menon 2009), where 7rp]!( s (e) is the exter¬ 
nal ring photon density per unit volume within the co¬ 
moving frame of the blob, 7/' is the angle between the 
scattering electron in the blob and the seed photon from 
the ring in the co-moving frame of the blob, dcr/de s dP s 
is the differential scattering cross-section, and dP and 
df2 e are the differential solid angles in the photon and 
electron directions, respectively. 

In order to evaluate 7i P h S ( e ) in the co-moving frame 
of the blob, we follow the method outlined by |Dermer fe 
Schlickeiser (1993), and make use of the following Lorentz 
invariants: n p h(e)/e 2 = 7T p h(e*)/e* 2 (where * denotes 
the host galaxy/stationary ring frame and unprimed no¬ 
tation refers to the co-moving frame of the blob) and 
e*/e = T(1 + /3/r), where T is the bulk Lorentz factor of 
the blob, /3 = Tibiob/c = \/l — F -2 is the blob’s velocity 
in units of c, and /i is the cosine of the angle that the 
incoming ring photons make with respect to the blob in 
the co-moving frame of the blob. The external ring pho¬ 
ton density (in photons cm^ 3 e _1 ) evaluated in the host 
galaxy frame can be shown to be: 


/ *\ A pnv , 
”' t(f > = Trf 


P h(e*) 


2 tt V 1 - PV* 


(Dermer & Schlickeiser 1993), where d = y/r 2 + z 2 and 
yU* = cos(0) = z/\/r 2 + z; 2 , and where 2 is the dis¬ 
tance of the blob from the center of the ring (see Fig¬ 
ure]^. Here, -/V p h(e*) is the photon production rate (in 


photons s 


-1 „-i 


) of the ring due to synchrotron radia¬ 


tion. If one further assumes that the ring is optically 
thin; we can write 


iV(e*) = 


4nr h 0 h syn (e*) dr 


where h 0 is the thickness of the ring in the 2 direction, 
and where r m i n and r max denote the inner and outer radii 
of the ring, respectively (see Figure [2]). Here, 7z syn (e*) is 
the synchrotron photon production rate of the ring (see 
Equation ([I])) evaluated in the host galaxy frame. It then 
follows that rip“ s (e) can be evaluated as follows: 


ring/ N _ _ 1 

P h 1 ; 2 ttc r 2 (l + / 3 / 7) 2 


x S 



l-Pn*) 



h syn (e*) dr 


( 5 ) 


A head-on approximation is adopted, in which the scat¬ 
tered photon direction is set equal to the electron direc¬ 
tion in the blob frame, which then results in the 
ing angle: cosi/j = + \J\ — n 2 ^J\ — /j , 2 cos 4> 


scatte r- 

|Dermer 
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TABLE 1 

Model Parameters 


Global Parameter Value 

7min 

2.0 x 10 3 

7max 

1.0 x 10 4 

z 

0.361 

^obs 

1.4° 

Baseline Flux optica i 1.5 (mjy) 

Baseline Flux - 7 — ray 6-5 (10 7 phot cm -2 s — x ) 

Blob Parameter 

Value 

T’blob 

0.01 (pc) 

^initial 

4 

Tfinal 

19 

•^initial 

- 0.1 (pc) 

•^final 

0.1 (pc) 

•®blob 

4.0 

^blob 

0.03 (G) 

P ■ 

-Mnj 

1.65 X 10 43 (ergs s _1 ) 

Ring Parameter 

Value 

^ring 

0.0 (pc) 

T’ min 

0.09 (pc) 

T’max 

0.18 (pc) 

h 0 

0.018 (pc) 

•Sring 

4.0 

-^ring 

0.12 (G) 

& Menon 2009 

), where /z s = (/i 0 bs - /3)/(l - P/iobs), and 

F'obs — C0S($ o bs 

), where 0 o b s is the angle between the jet 

axis and our line of sight (see Figure |2|) , and (j> is the 

azimuthal angle about the jet axis. In order to keep the 

numerical integration to a minimum, a delta-function ap- 

proximation of the differential scattering cross-section in 

the Thomson regime is adopted: 


dcr 


de s dfh 


= a T S(Ct s - Ct e )S(e s -^6)11(1-e) ( 6 ) 


(Dermer & Menon 2009), where fl s refers to the solid 
angle subtended by the scattered photons, e = ye(l — 
cos?/>), and H( 1 — e) is a Heaviside function given by: 



Fig. 3. — SED produced by the model during the orphan 7 -ray 
flare. The SED is composed of synchrotron (blue), synchrotron 
self-Compton (green), and external-Compton (red) emission com¬ 
ponents. Multi-wavelength data obtained by the SMARTS, Swift , 
and Fermi telescopes during the epoch of the orphan 7 -ray flare 
presented in Figure |T| are plotted in black for comparison. 


range of scattered photon and electron energies modeled 
in our code (see Table [l]), and we find: 0.003 < e < 0.05. 
This range of e falls well within the Thomson regime. In 
the future, the code will incorporate the effects of the 
reduced Klein-Nishina scattering cross-section at higher 
energies. Since j>dfl = / dfi f d c/>, Equations Q-® can 
be combined to yield a simplified version of the external 
Compton photon production rate: 


dec; (7 ) 


f gT^o 
( 2 t r 



dy 


1 

r 2 (i + Pn) 2 



«e( 7 ) 


x n syn (e*) (1 - cos-i/j) H(1 - e) , ( 8 ) 

where now e = e s /y 2 (l — cos ip) and fi = (fi* — /3)/(I — 

The distribution of electrons in the blob is evolved each 
time step by solving the Fokker-Planck equation: 


dn e { 7 ) 
dt 


= —b m e c 2 


d_ 

£>7 


( 7V (7)) + n 2 


( 9 ) 


Hi 1 - e) = 


0 for (1 - e) < 1 
| for (1 — e) = 1 
1 for (1 - e) > 1 . 


(Kardashe'dl962), where b is the electron radiative cool- 


( 7 ) 


As discussed in Dermer & Schlickeiser ([T993|, Equa¬ 
tion ([6]) is a valid approximation of the scattering cross- 
section provided e < 1. The delta-function S(e s — ye) in 
Equation ([ 6 ]) implies e = e s /y. The dimensionless scat¬ 
tered photon energy (e s ) is defined as e s = hv s /m e c 2 . 
The orphan flare (presented in Figure [Tj) we model with 
our code occurs at an observational frequency of t'obs ~ 
10 23 Hz. Due to the relativistic motion of the blob, 
the scattered photon frequency ( v s ) in the co-moving 
frame of the blob is related to the observed frequency by 
v a = r'obs /D, where D is a relativistic Doppler boosting 
factor given by D = 1/T(1 —/3 cos0 o bs)- With these quan¬ 
tities defined, we compute the limiting values of e for the 


ing rate in ergs ^ 1 s _i , and q is the electron injection 
rate per unit dimensionless energy in cm ~ 3 s” 1 . In gen¬ 
eral, b = 6 syn + b ssc + b ec \ however, here we assume that 
b ~ b syn + b ec because ssc is only a minor component in 
the SEDs produced by the runs of our code presented 
here (see Figure wl). In the future we intend to carry out 
a more sophisticated calculation that includes a proper 
treatment o f the ssc coo li ng, w hich, as pointed out by 
Zacharias & Schlickeiser (20121, can strongly influence 
the evolution of n e (y) in time-dependent calculations of 
the non-thermal emission from some blazars. For sim¬ 
plicity, we neglect aging in the ring electrons. The syn¬ 
chrotron cooling rate for the blob electrons is given by: 


^■=3 OTT (taSf?) (1 °> 

(Dermer fe Menon|2009 ). The external Compton cooling 
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rate is given by: 


b e c = 7, ccr T 


m e c z 


ring 

V 


(e) de ( 11 ) 


(Dermer & Menon 2009), where the upper limit of the in¬ 


tegration ensures that the cooling rate is evaluated within 
the Thomson regime. 

The injection rate q is computed by defining an injec¬ 
tion power: 


r 'Tmax / ^ 

p inj = j {'im e c 2 )(q r y~ s ) ^-7rrg lob ) d'y 


( 12 ) 


->• q = 


(2 - s)P i 


ln J 


(m e c 2 ) (|7rr3 lob ) ( 7 Lx - 7min) 


— (for s > 2 ). 


The injection power parameterizes a physical mechanism 
at work within the blob, for example, turbulence or dif¬ 
fusive shock acceleration, that continuously rejuvenates 
the aging blob electron distribution. This injection power 
is a free parameter within our model and has been set at 


Pw 


10 43 

Equation (|9j) has the following analytical solution: 


ergs s 1 in the calculations presented here. 


3 (7 >*) = < 


n Q 7 


for £ = 0 


n o 7 - s (l -£) s - 2 

'"''Sy [1 - (1 - 0 s ' 1 ] for £ < 1 


97 _ 

m e c 2 b(s— 1) 

-OH-i) 


97 


m e c 2 b(s— 1) 


for £ > 1 , 


(13) 


where £ = bt^m e c 2 (Kardashev 1962). The initial rela¬ 
tivistic electron number density, n 0 , of the blob is deter¬ 
mined by assuming equipartition between the magnetic 
and electron energy densities within the blob: 


o 2 

^blob 

87 r 


1 1 max. 

/ ( 7 TO e c 2 )(n 0 7 _s )d 7 

J 7min 

m 


2 -s 
/max 


■7 2 - 

1 min 


~zt (for s > 2). (14) 


After the photon production rates have been computed 
for a range of photon energies at each time step, the flux 
is computed in the observer’s frame by multiplying the 
sum of the production rates (htot = h syn + n ssc + n ec ) by 
the volume of the blob, the energy of the photons, and a 
(de/di') conversion factor: 


1 


D 


vF{v) ohs = —p X x ntot(e) x -vrrg lob 


M 

m e c 2 


(15) 


where Z is the redshift, c?l is the luminosity distance, 
and D is a relativistic Doppler boosting factor, given by 
D = 1/T(1 — /? cos0 obs ). 


3. NUMERICAL SIMULATION 

The numerical calculations were carried out using the 
BLAZE code (people.bu.edu/nmacdona/Research.html). 
A spherical blob of radius 7'bi 0 b = 0.01 pc accelerates 
from a bulk Lorentz factor of Tinitial = 4 to Tfi na i = 19, 
after which it continues to propagate at constant speed 
along the spine of the jet. The parameters governing the 
blob acceleration have been chosen to reflect the acceler¬ 
ation derived for the emission knot associated with the 
orphan 7 -ray flare in PKS 1510—089 (see Figure]!]) based 
on the changing rat e of rotation of the observea optical 
polarization vector (Marscher et al. 2008[ 2010). The 
blob begins its acceleration toward tne ring (located at 
-Aing = 0 pc in the simulation) at Zintiai ~ — 0.1 pc and 
reaches Tg na i at Zfi na i ~ 0.1 pc after passing through 
the ring. The blob is tracked out to a distance of 
~end ~ 0.4 pc, after which the simulation is terminated. 
Although in actuality the sheath (and by extension the 
ring) will have some relatively low velocity along the jet, 
for the purposes of the numerical simulation, we treat 
the ring as stationary. 

The dimensions of the ring used in the simulation, 
namely, its minimum (r m i n ) and maximum (r max ) radius 
and its thickness (h 0 ) (see Figure [ 2 ]), are listed in Table 
[I] The size of the sheath (and by extension the ring) is 
thought to be a substantial fraction of the cross-sectional 
width of the jet (|Zavala & Taylor|2004| |Walg et al.|2013 1 
and the values chosen tor r m ; n and r max reflect the width 
of the jet in PKS 1510—089. Since we are positing that 
the ring represents a shocked segment of a continuous 
jet sheath (see Figure [9] below), the thickness of the ring 
( h 0 ) in our model reflects the length scale over which par¬ 
ticle acceleration acts on a small section of the sheath, 
perhaps as it propagated through a standing shock up¬ 
stream of its current location within the jet. While this 
parameter is the least constrained in our model, h 0 sim¬ 
ply scales the external ring photon density up or down, 
d epending on our choice of its value (see Equation (|5b). 

Nalewajko et al. (2014) make a case against the ability 
ol a sheath ol plasma surrounding the jet to provide suf¬ 
ficient numbers of seed photons to produce the orphan 
flares seen in PKS 1510—089 (i.e., Figure [1} without the 
sheath itself outshining the observed flux. We find no 
conflict, however, between our model and observations 
of this blazar. As discussed in their paper, the observed 
luminosity of the sheath can be shown to be (their Equa¬ 
tion (C3)): 


i s h,obs - 2.7 x 10 47 erg s 1 

x r 4 h (r s iA h ) 2 , 


lOpc 


2 /p \ -2 
20 


(16) 


where r is the distance along the jet, T is the blob’s bulk 
Lorentz factor, and r s h and 0 s h are the bulk Lorentz fac- 
tor and opening angle of the sheath, respectively. As 
|Nalewajko et al.| point out, the strong dependence of 
L s h j0 bs on l’ s h means that for even mildly relativistic 
sheaths the observed luminosity can become quite large. 
The ring of sheath plasma in our model, however, is 
located upstream of the radio core in PKS 1510—089. 
In particular, the radio knot correlated with the orphan 
flare w as traveling down the jet at ~ 0.3 pc day -1 in our 
frame (Marscher et al. 2010). From Figure]!] it is appar- 
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ent that the onset of the orphan flare preceded the large 
optical flare associated with the knot passing the radio 
core by about 20 days. Assuming that the knot moved 
at this observed rate along the jet for those 20 days, that 
would place our ring 6 pc upstream of the radio core. If 
we further assume that the radio core is at roughly a scale 
of ~ 10 pc from the central engine, that would place the 
ring at roughly a distance of r ~ 4 pc along the jet. This 
locates our emission region between the two commonly 
considered regions of 7 -ray production in blazar jets, 
namely, the near-dissipation region (r ~ 1 pc), where the 
external photon fields are dominated by emission from 
a dusty torus and the BLR, and the far-dissipation re¬ 
gion (r > 10 pc) near the core seen at 43 GHz. At this 
distance, utilizing the dimensions of our ring as corre¬ 
sponding to the jet sheath, we compute a sheath opening 
angle of 6 > sh ~ 1.29° (where 9 sh = 6 > max - 6 > min ). Here, 
0 m i n and 0 max refer to the inner and outer angles sub¬ 
tended by the ring relative to the jet apex, which we 
take to be 4 pc upstream of the ring-see Figure [9] -and 
can be shown to be: 0 min = tan - 1 (r m i n /4 pc) ~ 1.29° 
and 0 max = tan - 1 (r max /4 pc) ~ 2.58°, respectively. We 
point out that this value for the sheath opening an¬ 
gle is in rough agreement with the sheath opening an¬ 
gle of 0 s h ~ 4.6° used by |D’Arcangelo et al.| (|2009|) to 
explain variability in the polarization ot the blazar OJ 
287. If we further assume that the sheath at this scale 
is only mildly relativistic (r s h ~ 2 ) and take our final 
value for the bulk Lorentz factor of the blob (r ~ 19), 
Equation (16) yields L s h,obs ~ 1.5 x 10 45 erg s -1 , which 
is two orders of magnitude smaller than the observed 
7 -ray luminosity for PKS 1510—089 during this flaring 
epoch (L 7 ~ 5.4 x 10 47 erg s -1 ) and is roughly a factor 
of 3.6 smaller than the observed synchrotron luminosity 


that amplifies the field strength through shock compres¬ 
sion (see Figure [9] below). We also approximate that the 
magnetic fields within the blob and the ring do not vary 
appreciably on the timescale required for the blob to pass 
through the ring (~ 100 days in the AGN frame). This 
requires that the opening a ngle of the jet < 10 °, which 
agrees with 9 ~ 1° found by fjorstad et al.| (j2005|). 

The numerical integrations were performed using linear 
grids of 15 radial and azimuthal zones and logarithmic 
grids of 60 intervals in e and in 7 . In order to convert the 
resultant flux into the observer’s frame, the j et is t aken 
to lie at an inclination angle of 9 0 ^ s = 1. 4° j jorstad et 


(L t /L s 


100 


5.4 x 10 45 erg s -1 ). 


The blob and ring electron populations are both mod¬ 
eled as identical power-law distributions, n e ( 7 ), with 
identical power-law indices s (listed in Table [Tb, in rough 
agreement with the observed 7 -ray photon in dex of r_= 
2.48 (a = T — 1 ; s = 2a + 1 ) reported in|Marscher et alT 


al. 20051), with a redshift of Z = 0.361 (Thompson et 
aE ISSfJ), and a luminosity distance of di, = 1913 Mpc. 
These values are based on the observed properties of PKS 
1510—089 and yield a projected scale of 5.0 pc mas -1 . To 
better fit the observed light curves presented in Figure [T] 
we include a baseline flux in each band, which is added 
to the variable emission produced by the model. These 
baseline fluxes, listed in Table [l] reflect the more slowly 
varying emission from the jet plasma and are approxi¬ 
mated to be constant. Finally, the time step is chosen 
so that 24 samplings of the SED are made by the end of 
the simulation of the blob passing through the ring. 

4. RESULTS 

Figure [3] shows an SED produced by the model during 
the peak of the simulated orphan 7 -ray flare. The three 
non-thermal emission components included in the cal¬ 
culation are differentiated by color. The model broadly 
reproduces the observed optical and 7 -ray peaks of PKS 
1510—089’s SED (shown in black), but does not repro¬ 
duce the observed optical-ultraviolet excess (known as 
the big blue bump), nor the observed X-ray flux levels. 
The former spectral feature is believed to be associated 
with emission from the accretion disk surrounding the 
central engine, while the latter may be ssc emission from 
regions farther down the jet (as indicated by the weak 
cor relation between the X-ray and 7 -ray light curves; 


see 


Marscher et al. 2010). Neither of these components 


are making the assumption of an unbroken power-law of 
electron energies in both regions of our model. The range 
of electron energies over which we assume the power-law 
applies ( 7 mi n through 7 max) is listed in Table H] 
Constant magnetic fields are taken to permeate both 
the blob (.Bbiob) and the ring (B r i n g) with the values 
used in the simulation listed in Table [T] As discussed 


Nishikawa et al. (20131, the jet sheath is thought to 


(2010). By setting the power-law index to this value, we are included in the calculation, hence these spectral fea¬ 
tures are absent from the synthetic SEDs produced by 
the model. 

The resultant 7 -ray and optical light curves produced 
by the blob’s passage through the ring are shown in Fig- 
urej4|( highli ghted i n red) overlaid upon the observed data 
froinjMarscher et al. (2010). In particular, the baseline 
flux level (as discussed in §3) is demarcated by the black 
dashed line in each plot, whereas the variable emission 
produced by the model is indicated by the black dotted 
line. The superposition of these two lines in each plot 
creates the synthetic light-curves shown in red. The syn¬ 
thetic light curves reproduce, in general, the observed 
orphan 7 -ray flare without exceeding the observed op¬ 
tical flux during the flare. The lack of a strong optical 
flare demonstrates the model’s key ability to produce or¬ 
phan 7 -ray flares. The peak flux in the simulated 7 -ray 
flare occurs when the blob is at a distance of z = 0.08 pc 
downstream of the ring. This is a result of the interplay 
between the Doppler boost gained by the relativistic mo¬ 
tion of the blob, the external ring photon density within 
the co-moving frame of the blob, and the effects of ra¬ 
diative cooling. It should be pointed out that there is a 
small optical “bump” in our synthetic light curve (lower 


have a higher electron density than the spine. Since we 
assume equipartition, this implies that B r i n g > Bbiob- 
In addition, a lower magnetic field within the blob is 
needed to keep the radiative losses experienced by the 
blob electrons out of the fast cooling regime (see Equa¬ 
tion (fl3|)); however, the value used in our simulation 
(Bbiob = 0.03 G) is at odds with the observed value of 


B 


jet 


0.1 — 1 G for PKS 1510—089 reported in Marscher 


et al. ( |2010 ). This difference in field strength might be 


reconciled by the fact that, as discussed above, we treat 
the ring as being located upstream of the radio core in 
PKS 1510—089. It is consistent for the magnetic fields in 
the jet upstream of the core to be smaller if indeed the 
core at 43 GHz is associated with a recollimation shock 
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days from orphan flare peak 




Fig. 5. — Comparative light curves highlighting the effect differ¬ 
ent motions of the blob along the spine have on the orphan flares 
produced by this model. In all cases, the flux for each flare has been 
normalized to the flux of that flare’s peak. For comparison, the fit 
to the orphan flare shown in Figure HI is reproduced in red. The 
blob’s initial (T initial) and final (Ffi na i)bulk Lorentz factor for each 
run are listed in the upper right portion of the plot (color-coded 
to the corresponding flare). It is apparent from the blue curve 
that if the blob decelerates as it approaches the ring, the resultant 
flare has a steeper onset time in comparison with its decay time. 
A similar effect (albeit with a distinct profile) is produced if the 
blob moves through the ring at constant speed (green curve). In 
contrast to the relatively symmetric flare profile (shown in red), if 
the blob approaches the ring with a greater velocity (black curve), 
the resultant flare has a more gradual onset. 

5. POLARIMETRIC SIGNATURE OF A JET SHEATH 


Fig. 4.- PKS 1510—089 7 -ray (upper panel) and optical (lower 
panel) fight curves from 2009 (black circles)—a smaller segment of 
Figure [TT The model hght curves are overlaid (in red) and consist 
of the superposition of the model baseline level of flux in each band 
(the dashed black line) and the variable emission produced by the 
model (the dotted black line). As with Figurell] the vertical arrow 
in the lower panel marks the time when a superluminal knot passed 
through the 43 GHz core of PKS 1510—089, with the horizontal bar 
representing the uncertainty in this time. 

panel of Figure El). This gradual increase in emission is 
a result of the Doppler boost gained by the blob’s accel¬ 
eration along the spine of the jet. 

The model’s ability to produce 7 -ray flares with steep 
rise and decay times (as illustrated in the upper panel of 
Figure |4]) reflects the anisotropic source of seed photons 
used in the simulation (i.e., the ring). This is in con¬ 
trast to other structured jet models that invoke isotropic 
seed photon fields from a continuous jet sheath. The 
shape of the orphan 7 -ray flare shown in Figure [4] is very 
sensitive to the nature of the blob’s acceleration toward 
and through the ring. To further explore the effect the 
blob’s motion along the jet has on the orphan 7 -ray flares 
produced by this model, the parameters governing the 
blob’s acceleration along the spine of the jet were varied, 
namely, the blob’s initial (Tinitiai) and final (F fi na i ) bulk 
Lorentz factors. In Figure[5j a series of comparative light 
curves is shown illustrating the effect that different blob 
accelerations (or lack of acceleration) have on the orphan 
flares produced by this model. It is apparent from Fig¬ 
ure [5] that by varying the blob’s acceleration through the 
ring, this model is able to produce orphan flares both 
with steeper and more gradual onsets and decays. The 
differences in these light curves highlight the potential to 
use the observed shape of orphan 7 -ray flares to infer the 
nature of the acceleration of the emission region through 
blazar jets upstream of the 43 GHz core. 


The presence of a sheath of plasma surrounding the 
spine of a blazar jet is a central component to our model 
of orphan 7 -ray flares. Several authors have discussed 
what the observational signature of a jet shea th within 
a blazar would be (see Pushkarev et al. 2005 and refer¬ 
ences therein). They posit that a jet sheath represents a 
shear layer between the plasma of the jet and the ambi¬ 
ent medium into which the jet propagates. In this sce¬ 
nario, the initially helical or chaotic magnetic field lines 
of the outer layers of the jet are “stretched” out due 
to the velocity gradient between the jet and the ambi¬ 
ent medium. This shear, if it exists, would result in a 
buildup of magnetic field on the outer edges of the jet, 
with the orientation ofjthat field aligned nearly parallel 
to the jet boundary (Wardle et al. 1994). The direction 
of the magnetic field within a blazar is inferred from the 
orientation of the EVPAs measured by the VLBA, which, 
barring t he effects of opacity and relati vistic aberration 
(see Lyutik ov, Pariev, fe Gabuzda|2005 ), are assumed to 
be perpendicular to the projection of the magnetic field 
onto the plane of the sky. Therefore, the observational 
signature of a jet sheath would be an increase in the frac¬ 
tional polarization toward the edges of the jet, with the 
EVPA perpendicular to the jet boundary, and therefore 
roughly perpendicular to the jet axis. We point out that 


an altern ative interpretation (as discussed in Lyutikov. 


Pariev, & Gabuzda 2005) is that the EVPA being" per¬ 
pendicular to the jet boundary is the result of the jet 
flow carrying a large scale helical magnetic field. A mea¬ 
sure of t h e Faraday rotation gr adient , if it exists, (see 


Gabuzda, Reichstein, & O’Neill 
of thejet of PKS 1510—089, wou 


2013) across the width 
d potentially be able to 


distinguish between these two scenarios. We adopt the 
prior interpretation for this paper. 

A complication in detecting the polarimetric signature 
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Fig. 6. — A stacked map of 43 GHz images of PKS 1510—089 spanning twenty epochs of observation from 2008 to 2013. The black 
contours correspond to total intensity (contour levels are (mjy beam -1 ): 2.3 4.6 9.2 18.4 36.8 73.6 147.2 294.4 588.8 1177.6), whereas the 
underlying color scheme corresponds to polarized intensity (see color bar to the right for the flux levels) with the EVPAs denoted by black 
line segments. All images have been convolved with a Gaussian beam, shown on the lower right of the stacked map. An arbitrary jet 
spine has been plotted in red, across which ten transverse slices through the data are taken (shown in blue). The profiles of the emission 
parameters along the first slice are shown in Figure [7] The shaded gray regions in the above map highlight the nominal location of the jet 
sheath of PKS 1510—089 as determined by the procedure discussed in §5 and highlighted in Figure [8] 


of a jet shea t h wit hin a blazar is discussed by IClausen/ 
Brown et al. (2013). Due to the effects of Doppler beam¬ 
ing, the relativistic spine of the jet (r sp j ne ~ 20) is intrin¬ 
sically far brighter than the surrounding sheath (which 
is usually assumed to be non- to mildly-relativistic: 
r s heath ~ 2). To overcome this vast difference in intrinsic 
brightness, several authors have imple mented a method 


of “stacking” radio map s (see, e.g., Fromm et al. 2013, 
Zamaninasab et al.|2013). This stacking method involves 


aligning individual radio images of a blazar based on the 
location of the radio core (which is assumed to be an es¬ 
sentially stationary feature within the jet at a given fre¬ 
quency). The images are then added together, with the 
expectation that noise and transient phenomena (such 
as superluminal blobs) will smooth out in the stacked 
map, and, in contrast, any standing features within the 
jet (such as a jet sheath) will be highlighted as more im¬ 
ages are added to the stacked map. With this in mind, 


we stacked VLBA total intensity and polarization images 
of PKS 1510—089 at 43 GHz in order to see whether the 
polarimetric signature of a jet sheath is present in this 
blazar. In particular, we selected images from twenty 
epochs of observation from 2008 to 2013 made with 
VLBA data from the VLBA-BU-BLAZAR pro gram (see 
www.bu.edu/blazars/VLBAproject.html and Jorstad et 


al. (2013) for a general description of the observations 
and data analysis). The chosen epochs correspond to 
periods of relatively quiescent jet activity. Our stacked 
map for PKS 1510—089 is shown in Figure [6] We do see 
an increase in the degree of polarization toward the edges 
of the jet. The EVPAs plotted in Figure [6] are inclined 
to the jet spine axis (highlighted in red) and are predom¬ 
inantly perpendicular to that axis (as predicted by the 
sheath scenario) in both the southern half and the up¬ 
per northern half of the jet. We examine ten transverse 
slices through the data (shown in blue). The variations 
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Fig. 7.— Variations of the emission parameters transverse to 
the jet axis (dashed red line) for slice 1 (shown in Figure 6): to¬ 
tal intensity (top panel), polarized intensity (upper middle panel), 
fractional polarization (lower middle panel), and EVPA (bottom 
panel). The fractional polarization increasing toward the edges of 
the jet is a predicted polarimetric signature of a jet sheath (as 
discussed in §5), with the EVPA transverse to the axis. 

of the Stokes parameters along the first slice are shown 
in Figure [ 7 ] We indeed find that the fractional polariza¬ 
tion increases toward the edges of the jet, as predicted 
by the sheath scenario. These two features, namely, the 
increase in fractional polarization and the orientation of 
that polarization perpendicular to the jet axis, are con¬ 
sistent with a sheath in PKS 1510—089. The variations 
of the Stokes parameters along the other nine slices are 
qualitatively similar to the first, with the exception that 
after the fifth slice we no longer see signs of a sheath 
component in the southern half of the jet. 

The total and polarized intensity profiles along the 
slices (highlighted in Figure [7| are used to determine 
the nominal location of the sheath within the jet, as il¬ 
lustrated in Figure [H] In particular, we align the total 



Distance from Jet Spine [mas] 


Fig. 8.— Illustration of the procedure used to determine both 
the location of the sheath within the jet and its contribution to 
the observed emission by utilizing the total and polarized intensity 
profiles (shown in FigureUll of the slices taken through our stacked 
map (shown in Figure fbjL The polarized intensity profile along 
slice 1 (dashed blue line corresponding to the right-hand axis) is 
overlaid upon the total intensity profile for slice 1 (solid blue line 
corresponding to the left-hand axis). We use this double-peaked 
profile to highlight the likely location of the sheath within each slice 
(shaded gray regions). The sheath’s contribution to the jet’s total 
intensity profile is estimated by the values of total intensity (red 
circles) that are cospatial with the peaks of the polarized intensity 
profile (demarcated by dashed vertical lines). This procedure is 
then repeated for each slice, thus tracing out the sheath (shaded 
gray regions) shown in Figure |6| 

intensity and polarized intensity profiles for each slice 
and use the double-peaked nature of the polarized in¬ 
tensity (a signature of the jet sheath) to highlight the 
location of the sheath within each slice (the shaded gray 
regions shown in Figure [8]). The width of the sheath on 
each side of the slice is cletermined, arbitrarily, by the 
location where the polarized flux falls below 0.85 times 
the peak value (the solid black vertical lines in Figure 
[8]) . This percentage yields a conservative estimate of the 
extent of the sheath. The total intensity of the sheath 
plasma on each side of the jet is approximated by the val¬ 
ues of the total intensity profile (the red circles in Figure 
[8j) that are aligned with the peaks of the polarized inten¬ 
sity profile (the dashed black vertical lines in Figure [8]). 
We then repeat this procedure for each slice throughout 
our stacked map, thus tracing out the sheath (shaded 
gray regions shown in Figure pi). 

An estimate of the bolometric luminosity (Tboi) of the 
sheath in the observer’s frame is computed by adding 
all of the flux contained within the gray shaded regions 
of our stacked map (Figure [6]). Based on the distance 
to PKS 1510—089, we convert this inferred sheath flux 
at 43 GHz into a spectral luminosity (L u ) and assume 
a power-law of the form L v = Cv~ a . After solving for 
the constant of proportionality (C ), we then integrate to 
obtain: 


Abol — 


L„ dv ~ 3 x 10 45 erg s 


-l 


(17) 


where we have assumed a spectral index of a ~ 1.0 
(in rough agreement with the mean spectral index of a 
lar ge sample of flat spec trum radio quasar jets found 
by Hovatta et al. 20141 and limits of integration of 
Fmin = 10 y Hz and iz max = 5 x 10 13 Hz. It is inter¬ 
esting to note the similarity between this value and the 
value obtained i n §3 for the upstream sheath luminosity 
using Equation (16). The above estimate indicates that 
the sheath within PKS 1510—089 is potentially a very 
important source of seed photons even at parsec scales 
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Fig. 9.— A schematic of the relative locations along the jet of both the ring of shocked sheath plasma in our model and the location 
of the radio core/sheath detected further downstream in the stacked radio images of PKS 1510—089 shown in Figure [6] As discussed in 
§3, we posit that the ring is located ~ 4 pc from the central engine while the radio core in the stacked map is located further downstream 
at a scale of > 10 pc from the central engine. We propose that the radio core in PKS 1510—089 is associated with a recollimation shock 
that compresses initially tangled magnetic field along the spine of the jet (shown in red) an orders that field perpendicular to the jet axis 
(the red vectors just to the right of the recollimation shock). In contrast to the spine, velocity shear between the sheath and the ambient 
medium (blue vectors denote relative speed) aligns the magnetic field lines on the outer edges of the jet to be roughly parallel to the jet 
axis resulting in the spine-sheath polarization signature we detect in our stacked map of PKS 1510—089 (shown to the right). 


along the jet. As a further comparison, we compute the 
bolometric luminosity (in the observer’s frame) of the 
ring of sheath plasma in our model assuming, in con¬ 
trast to the simulation, that the ring moves along the 
jet with a bulk Lorentz factor of r r i ng = 2. We find 
L™f ~ 3 x 10 44 erg s _1 , roughly an order of magnitude 
smaller than the above sheath estimate. 

Finally, a crude estimate of both the inner (r m ; n ) and 
outer (r max ) radii of the sheath was made by averaging 
the values obtained for these parameters relative to the 
jet spine (demarcated by the dashed red vertical line in 
Figure| 8 | for each slice through the stacked map. We find 
that r m i„ ~ 0.4 pc and r max ~ 0.8 pc. These values are 
roughly a factor of 4.5 larger than the inner and outer 
radii used in our model of the ring (see Table [l]) that 
we posit to exist upstream of the region highlignted in 
Figure [ 6 } It is interesting to note that the opening angle 
of the sheath used in our model (# s h = 9 max — 0 m i n ~ 
1.29°; 0 m i n = 1.29° and 0 max = 2.58° as discussed in §3) 
predicts that r m ; n = z tan(f? m j n ) ~ 0.4 pc and r max = 
z tan( 0 max ) ~ 0.8 pc for the sheath at a distance of 
0 = 17 pc from the central engine, which is roughly the 
scale where the radio core is situated in our stacked map 


at 43 GHz (Marscher et al. 


Although the plasma m t 


2010 ). 


lis sheath needs to move rel- 


ativistically for beaming to allow us to detect it, given 
this transverse velocity gradient, we infer the presence 
of an even slower sheath upstream, where the ring in 
our model could be located (as discussed in §3 and il¬ 
lustrated in Figure^). An enhancement, such as a ring 
of shocked plasma, m this more slowly moving sheath of 
plasma could potentially supply the seed photons that 
our model requires to explain the observed orphan flares 
within PKS 1510-089. 


6. SUMMARY AND CONCLUSION 

We have developed a model to explain the origins of 
orphan 7 -ray flares from blazars and have applied this 
model to a specific orphan flare within the blazar PKS 
1510—089. In this model, a shocked segment of a slow jet 
sheath provides a localized source of seed photons that 
are inverse-Compton scattered by electrons contained 
within a blob of plasma moving relativistically along the 
spine of the jet. This inverse-Compton scattering results 
in the production of a sharply peaked orphan 7 -ray flare. 
Synthetic light curves produced by this model are able 
to reproduce the general features of the observed 7 -ray 
light curve during an orphan flare event associated with 
the ejection of a superluminal knot from the base of the 
jet of PKS 1510—089. The model is also able to repro¬ 
duce the double-peaked profile seen in the SEDs of most 
blazars. This model can potentially be used to discern 
the nature of the motion of the blob along the spine of 
the jet based on the shape of the observed orphan flare. 
We have stacked radio images of PKS 1510—089 obtained 
with the VLBA at 43 GHz at multiple epochs. We find 
the polarimetric signature of a sheath in our stacked map, 
thus lending observational support to the plausibility of 
this model of orphan 7 -ray flares within blazars. In the 
future, we plan to extend our modeling efforts in an at¬ 
tempt to reproduce orphan 7 -ray flares observed in the 
blazars 3C 273, 3C 279, and 4C 71.07. We also plan to 
create stacked maps of the jets of these objects to in¬ 
vestigate whether we can detect additional spine-sheath 
polarimetric signatures. 
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